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ABSTRACT 


' —  A  procedure  for  deriving  the  variance  of  the  forecast  error  for  Winters’ 
Additive  Seasonal  Forecasting  system  is  given.  Both  point  and  cumulative 
T>step  ahead  forecasts  are  dealt  with.  Closed  form  expressions  are  given  in 
the  cases  when  the  model  is  (i)  trend-free  and  (ii)  non-seasonal .  The  effects 
of  renormalization  of  the  seasonal  factors  is  also  discussed.  The  fact  that 


the  error  variance  for  this  system  can  be  infinite  is  discussed  and  the  rela¬ 
tionship  of  this  property  with  the  stability  of  the  system  indicated.  Some 
recommendations  are  given  about  what  to  do  in  these  circumstances. 


Keywords:  Winters  Additive  Seasonal  Model;  Forecast  Error  Variance; 


Confidence  Intervals;  System  Stability 
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All  forecasts  are  wrong.  As  a  result  It  Is  rarely  adequate  to  give  a 
forecast  by  itself.  Some  measure  of  the  reliability  or  accuracy  of  the  fore¬ 
cast  is  also  required.  Such  a  measure  is  usually  provided  by  a  confidence 
interval.  This  is  an  interval  within  which  the  future  value  we  are  forecast¬ 
ing  will  fall  with  a  prescribed  probability.  The  philosophy  and  construction 
of  such  intervals  are  well  documented  in  a  number  of  forecasting  texts  and  we 
will  not  consider  them  further  here.  See,  for  example.  Brown  (1962), 

Montgomery  and  Johnson  (1976)  and  Bowerman  and  O'Connell  (1979). 

We  simply  note  that  an  essential  ingredient  in  all  the  expressions  and 
formulae  involved  in  such  intervals  is  the  variance  of  the  forecast  error. 
Consider  the  T-step  ahead  forecast  made  at  time  t  ,  i.e.  the  forecast  of 

A  ^ 

made  at  time  t  ,  X^(T)  say.  The  corresponding  error  is  e^(T)  *  X^^^  -  X^(T). 

2  2 

We  shall  denote  its  variance,  var[e^(T)],  by  Oj  .  It  is  clear  that  Oy  must 

A 

play  a  central  role  in  any  discussion  of  the  usefulness  of  X^(T)  in  predicting 
^t+T  * 

2 

The  value  of  Oj  can  usually  be  computed  on  the  assumption  that  the  under¬ 
lying  model  is  valid.  By  underlying  model  we  mean  the  model  which  the  forecast¬ 
ing  system  assumes  is  generating  the  data.  The  variance  can  certainly  be  derived 
for  the  ARMA  models  of  Box  and  Jenkins,  and  the  interested  reader  is  directed  to 
their  book  (1976).  It  can  also  be  derived  for  all  the  Exponential  Smoothing 
systems  discussed  by  Brown  (1962)  and  others  considered  more  recently  by 
Sweet  (1981). 

However,  no  results  are  available  for  Winters'  seasonal  systems.  This  is 
particularly  unfortunate  since  these  systems  are  amongst  the  most  commonly  used 
in  practice.  They  are  relatively  simple  to  implement  and  intuitively  appealing. 
Both  systems  are  described  in  detail  in  most  forecasting  texts,  e.g.  Montgomery 
and  Johnson  (1976),  Thomopoulos  (1980),  and  Bowerman  and  O'Connell  (1979).  In 


this  last  work,  approximate  confidence  intervals  are  given  for  Winters'  systems. 
As  we  shall  see,  however,  these  intervals  poorly  reflect  the  behaviour  of  the 
error  variance. 

The  Winters'  Multiplicative  seasonal  form  is  a  non-linear  system  and  it  is 
difficult  to  see  how  any  useful  information  about  the  forecast  error  can  be 
obtained  directly.  Approximation  and  simulation  appear  to  be  the  most  sensible 
tools  here.  In  the  case  of  the  additive  model,  however,  we  can  obtain  some 
results. 

The  purpose  of  this  paper  is  to  give  a  procedure  for  deriving  the  variance 

of  the  forecast  error  for  the  Winters'  Additive  Seasonal  forecasting  system.  We 

deal  with  not  only  point  forecasts  of  future  values,  ,  but  alsp  cumulative 

forecasts,  i.e.  forecasts  of  Y.  »  “  t  .  ,  made  at  time  t  . 

t,i  I 

Expressions  are  given  for  the  variances  of  the  individual  components  of  the 
forecasts,  i.e.  level,  trend  gradient  and  seasonal  factors,  and  for  the  covari¬ 
ances  between  them.  This  allows  the  construction  of  the  variance  of  any  par¬ 
ticular  forecast  and  hence  of  the  corresponding  forecast  error.  It  also  enables 
us  to  construct  confidence  intervals  for  the  components  themselves,  and  linear 
combinations  of  them. 

Two  special  cases  are  considered,  viz.  when  there  is  no  trend,  and  when 
there  is  no  seasonality.  This  latter  is  the  well-known  Holt-Winters  non-seasonal 
forecasting  system.  In  both  cases,  closed  expressions  are  obtained  for  the 
variances.  In  the  completely  general  seasonal  case,  the  expressions  derived 
involve  a  number  of  unknowns  which  are  obtained  by  solving  a  set  of  linear 
equations. 

Some  necessary  discussion  is  also  given  about  the  stability  of  these 
forecasting  systems. 
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1.  THE  RESULTS 

The  underlying  model  assumed  by  the  Winters'  additive  seasonal  system  1$ 


given  by 


»  m  +  bt  +  S|^  + 


where  t  =  rn  +  k  ,  and  n  Is  the  length  of  the  season.  The  values 
{S|^  :  k  »  l,2.....n}  are  the  additive  seasonal  factors  and  {a^}  Is  a  sequence 
of  Independent  Identically  distributed  random  variates  of  zero  mean  and 
variance  o'  . 

The  standard  form  of  the  corresponding  forecasting  system  Is  given  by 

'"t  “  ®to^^t‘^t-n^  *  ^ 

bt  =  oiCvVi)  ^  ^2) 

^t  *  ®2^^t“"*t^  ^^’“2^^t-n  • 

The  T-step  ahead  forecast  Is 

\(V  *  +  Tb^  ^  (3) 

where  T  *  rn  +  k,  (k  =  l,2,...,n;  r  ^  0)  . 

The  cumulative  forecast,  I.e,  the  forecast  of 


^t  T  ’  I  h*i  ^t  T  ■  i  *t^^)  * 

’  1=1  1=1  ^ 

Our  derivation  of  the  forecast  error  variances  parallels  the  development 
given  for  Exponential  Smoothing  models  by  Brown  (1962).  By  definition,  X^(T) 
depends  on  only  past  values  of  a^  and  so  Is  Independent  of  .  Thus, 

Oj  ■  o  +  var[Xj.(T)]  .  Further,  Xj(T)  can  be  expressed  as  a  linear  combination 
of  these  past  values  of  a,,  and  so  its  variance  Is  proportional  to  .  Thus, 


^  *  o^Cl  *  V(T)]  .  When  dee liny  with  cufflulative  forecasts,  the  same  aryument 
^h^4%  iTIie  error  now  is  j  ’  ^t.T  variance 

of  the  form  o^CT  +  Vg(T)]  . 

In  practice.  0^  is  unknown.  It  may  be  estimated  from  the  data,  however, 

9  2 

using  the  fact  that  o  [1  +  VU)]  *  Oj  .  the  variance  of  the  one  step  ahead 

forecast  error.,  ,  Thus,  we  can  estimate  Oj  directly  from  the  forecast  errors, 

2 

and  knowledge  of  V(l)  yields  an  estimate  of  o  . 

,  We  derive  ..V(T)  and  Vg(T)  from  the  variance-covariance  matrix  of  the 
estimates  of  the  individual  components  of  the  forecast  i.e,  m^,  and 

k  =  1.2»...,,n  .  These  variances  are  of  some  interest  in  their  own 
right  for  the  construction  of  confidence  intervals  for  the  components,  and  we 
display  them  below.  Details  of  the  derivations  are  given  in  Appendix  1. 

®m  *  ~  ^  (n-i)d^]  ♦  2(l-oj)(dy-d^)}o^ 

%  =  Var(b^)  .  2aJ<^(do-d„)o2 

®s  ’  “  2o|(l-<\))^(V‘^l)o^.  k  -  1.2 . n 

®mb  *  Cov(m^.b^)  =  <|2aj(2-aj)(dQ-d^)a^ 

®bk  ’  *  ^k-1  ’  ‘‘nrk+P®^*  "  '  . "  • 

®mk  ’  ‘^°''^'"t»®t+(c-n^  *  ‘»0«2^1-‘\))t(l-Oi)(d^_^-d^)  +  k  =  1,2 . n 

*  ®ji  “  ^®''^Vi-n»*t+j-n^  ■  '  ‘‘j-i-1  ‘  ‘^j-i+1^®^*  ^  ^  <  j 

Before  discussing  the  derivation  of  the  values  (dQ,d^,... ,d^}  we  note  the  form$ 
of  V(T)  '  iind  V-(T)  . 


Now,  V(T)a^  »  var(ra^  +  Tb^  +  ^t+k-n^ 

2  2  2  2 

‘  2(On,k  "■  '^%b  *  Hk^  ' 

Thus, 

V(T)  -  2oi§a^(dQ-d^)T2 

+  2aoaj[«^(2-Oj)(dg-d^)  +  agd-a^) (d^_^  -  d^  +  d^_j  -  ‘‘n-k+1^^^ 

+  *  (‘^k-rVk+l^^ 

+  2a^(l-ai)(dQ-d^)  +  2o|(l-e^)^(dQ-dj)  +  c^c^Lnd^  +  2  (n-i)d^]  (4) 

In  the  same  way,  using  the  fact  that 

h.T  T(T*l)bj  ♦  , 

we  can  show  that  when  T  =  pn  +  q,  p>,0,  l£q£n; 

♦  2'^(l-'>o)^C(l>*l)^(<lo  <lq)  +  P^fVVq*  *  P(P*‘)(VWq''‘o)^ 

*  «0«1»2(1-“0>I:<2P*»<V''-.>  *  (''n.q-'*q>]T(Tn) 

♦  Z«0O2(l-'>0)Cp(2-«l)(V''n>  *  <W‘'n*^n-q’ 

*  «1  J,  (“k-VklJT  (5) 

Now  the  numbers  dQ,d^,...,d^  are  necessary  to  evaluate  the  variances. 

These  {n*l)  values  are  derived  by  means  of  an  algorithm  which  Is  given  and 
discussed  In  Appendix  2.  It  Is  a  fairly  straightforward  algorithm  due  to 
Wilson  (1979).  It  has  the  property  that  It  also  checks  the  stability  of  the 
forecasting  system.  As  we  shall  see  In  Section  3  this  Is  most  useful  for  this 


forecasting  system.  Before  discussing  any  of  the  technical  aspects  of  the 
procedure  described  here,  we  present  some  particular  cases  which  have  closed- 


form  solutions. 


2.  SPECIAL  CASES 


2.1.  Trend-free  additive  seatonal  model. 

In  this  case  the  model  assumed  Is  a  constant  level  plus  addi¬ 
tive  seasonals.  The  forecasting  equations  are  obtained  by  removing  the  second 
equation  from  the  original  three  (2).  and  dropping  b^^^  from  the  first. 

Here  the  values  ^g*^p****<^n  ^  derived  algebraically.  The  variance- 

covariance  matrix  of  the  components  Is  now  obtained  from 

al  -  nagh  ^ 

2 

®s  "  ®2(l"«^)nn-l)<i^  ♦  » 

"ink  *  >  '  <  I"  1"  • 

•ij  ■  »j1  ■  -  ".k  •  l<1<J<n. 

2 

where  h  »  a  * 

Thus, 

.  ‘  -  O-^oXi-oz*  , 

r*  n-<\,')M'-oi2)  • 

It  Is  Interesting  to  note  that  for  this  model  V(T)  Is  Independent  of  T  . 

For  the  cumulative  forecast  error,  we  find  that  for  T  »  pn  +  q,  (p  ^  0, 

I  <1  <  n).  ♦  c|(l-c^)^(pT  ♦  pq  +  q) 

where  P  •  oiQCn<iQ  ♦  O2(l-oyj)]h/o^ 

2.2,  Non-soasonal  Holt-UInters*  system. 

The  appropriate  equations  are  obtained  by  removing  the  third  equation 
from  (2)  and  S^^^^  from  the  first.  This  Is  a  well-known  non-seasonal  fore¬ 
casting  system  and  often  used  for  linear  trend  models.  It  Is  a  direct  compet¬ 
itor  of  Brown's  second  order,  pr  double,  exponential  smoothing  which  we  will 
consider  shortly.  The  underlying  model  Is  assumed  to  be  a  simple  linear  trend 

1.0.  s  Is  removed  from  (1). 
k 
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We  obtain 


®m  •  3o|Qoi^)hj  , 

®fnb  "  ‘\)“l<2-<*i)»‘i  . 

2 

where  h^  =  o  /(4-2o^-o^Oj)  . 

We  can  deduce  In  the  usual  way  that 

V(T)  =*  [2agO^T^  +  2o^a^(2-a^)T  +  (2o^+2oj-3o^aj) ]hj/o^  (7) 

and 

Vg(T)  »  [ajjajT^  +  ■*■  (4o^+4aj-2oiQOj-oiQOj) ]T^hj/2o^ 

2.3.  Brown's  Double  Exponential  Smoothing. 

It  Is  well-known  that  this  system  is  equivalent  to  the  Holt-Winters' 

2 

non-seasonal  form  above  If  we  choose  >  1-6  and  (l-6)/(l'«'6)  where 

6  *  1-a  and  a  is  the  smoothing  constant.  Thus,  the  corresponding  results 

can  be  obtained  by  making  these  substitutions  In  the  equations  of  Section  2.2. 

2  2 

above.  The  results  for  o^,  o^jj  and  V(T)  may  be  found  In  Brown  (1962) 
and  Vg(T)  Is  given  In  Bowerman  and  O'Connell  (1979). 

2.4.  Continuously  re-normalized  seasonal  factors. 

In  the  general  additive  seasonal  model  given  by  (1)  it  is  usually 
n  • 

assumed  that  I  s,^  *  0  .  This  Is  done  to  ensure  some  measure  of  independence 
■k-»l  ■  ' 

between  the  level  of  the  data  and  the  seasonal  pattern.  It  Is  always  recom¬ 
mended  that  when  the  corresponding  forecasting  system,  as  given  by  (2),  Is 
1nH>1emented,  the  seasonal  factors  sum  to  zero  Initially.  What  Is  to  be  done 


8 


thereafter  Is  not  so  dear.  Because  of  the  revision  equation  for  the  seasonal 
factors  they  will  no  longer  sum  to  zero  after  the  first  observation.  We  can 
renormalize  the  seasonal  factors  at  any  time  by  subtracting  from  each  the 
average  of  the  set,  I.e.  the  most  recent  n  •  Whether  we  should  do  so  or  not 
Is  not  clear.  Some  authors  seem  to  recommend  against  It,  e.g.  Bowerman  and 
O'Connell  (1979);  some  regard  It  as  an  optional  modification,  e.g.  Montgomery 
and  Johnson  (1976);  some  recommend  renormalization  once  per  season,  e.g. 
Chatfleld  (1978);  and  some  recommend  continuous  renormalization,  I.e.  after 
every  revision  of  seasonal  factors,  e.g.  Thomopoulos  (1980). 

Our  only  Interest  here  Is  In  how  such  a  procedure  affects  the  forecast 
error  variance.  In  the  cases  of  occasional  or  purely  seasonal  renormalization 
the  situation  Is  very  complex  and  we  have  nothing  to  say,  except  that  the 
effect  on  the  error  variance  appears  to  be  small.  In  the  case  of  continuous 
renormalization,  however,  the  following  result  may  be  applied.  An  outline  of 
the  proof  Is  given  In  Appendix  1,  Continuous  renormalization  of  the  general 
forecasting  system  given  by  (2)  yields  exactly  the  same  forecasts  (and  so 
forecast  errors)  as  running  the  system  without  renormalization  but  replacing 
oio»<»i»«2  «5»®i  respectively,  where  aj  »  -  «»2(l-o^)/n  , 

oj  »  ,  and  aj  ■  a2/(l+a2/n)  .  As  a  consequence,  the  appropriate 

covariance  matrix  elements  and  values  of  V(T)  and  V^(T)  can  be  obtained 
as  In  Section  1  above  by  replacing  o^io^toi^  8y  oiQ,mpa^  respectively.  The 
same  holds  true  for  the  trend-free  version  discussed  In  Section  2.1. 


3.  TECHNICAL  CONSIDERATIONS 
3.1.  Infinite  Error  Variance  and  Stability. 

In  smoothing  systems  in  general  the  variance  of  the  forecast  error 
increases  as  the  values  of  the  smoothing  constants  increase.  For  example,  in 
the  case  of  the  simple  exponentially  weighted  moving  average  (SEWMA), 

A  A 

=  aX^  +  (l-a)X^_j  ,  the  variance  of  the  forecast  error  is  proportional  to 
2/(2-a)  which  clearly  increases  with  o  in  (0,1).  The  need  for  higher  values 
of  o  in  practice  reflects  the  fact  that  the  underlying  level  is  changing 
rapidly.  Consequently,  a  more  responsive  forecast  is  needed.  Moreover,  the 
inherent  instability  of  the  underlying  model  is  reflected  in  the  increase  in 
forecast  error  variaoce.  This  variance  is  finite  while  a  remains  in  the 
stability  region  of  the  system,  i.e.  (0,2)  for  the  SEWMA.  Note  that  the  set 
of  values  from  which  o  is  usually  chosen  is  a  subset  of  the  stability  region 
of  the  system.  Thus,’  a  stable  system  always  results  for  the  SEWMA,  This  is 
also  true  for  all  the  General  Exponential  Smoothing  models.  However,  it  is 
not  true  for  the  seasonal  system  under  consideration  here. 

It  is  a  somewhat  surprising  and  problematic  fact  that  there  are  choices 
of  the  smoothing  constants  Oq,  o^,  012  usual  range,  (0,1),  which 

yield  an  infinite  variance  for  the  forecast  error.  The  algorithm  given  in 
Appendix  2  checics  for  this  possibility  which  indicates  that  the  forecasting 
system  is  unstable.  If  such  a  situation  arises  clearly  no  meaningful  confi¬ 
dence  intervals  can  be  constructed.  More  importantly,  however,  we  must  decide 
how  to  interpret  this  knowledge  of  the  system's  instability.  The  concept  of 
stability  for  a  systan  of  difference  equations  is  an  important  one  but  rarely 
discussed  in  the  context  of  forecasting  systems.  Two  useful  exceptions  to 
this  are  the  papers  by  McClain  (1974)  and  Brenner  et  al,  (1968). 


In  essence,  a  forecasting  system- SliCh  as  (2)  above  Is  stable  if  the 
influence  of  earlier  observations  decreases  with  the  passage  of  time.  Thus, 
the  forecasts,  or’ (equivalently)  the  forecast  errors,  are  Influenced  more  by 
recent  observations  than  by  those  In  the  past.  It  Is  worth  noting  that  this 
Is  also  the  essential  philosophy  of  exponential  smoothing  systems,  and  Indeed, 
mostt  forecasting  systems. 

The  Converse  of  this  Is  that  In  an  unstable  system  past  observations  have 
a  constant  or  even  growing  Influence  on  future  forecasts.  As  an  Illustration, 

A 

consider  the  SEWMA  again.  The  forecast  can  be  written  as  a  weighted 

average  of  all  past  observations.  At  time  t  ,  the  weight  given  to  Is 

'k  • 

a(l-o)  .  Clearly,  if  o  lies  outside  the  stability  region  (0,2)  this  weight 
Increases  with  '  k  so  that  data  In  the  most  remote  part  have  greatest  Influ¬ 
ence  upon  the  forecast'.  Equally,  If  a  »  2,  all  observations,  however  distant 
In  time',  nlake  the  same  contribution  to  the  forecast. 

the  Inescapable  concluSI dh  Is  that  It  would  be  extremely  unwise  to  select 
smoothing  constants  which  do  not  Tie  In  the  stability  region  of  the  system. 

On  the  one  hand,  we  can  tolerate  a  high  (but  finite)  forecast  error  variance 
because  this  represents  a  trade-off  between  accuracy  and  robustness.  In  an 
effort  to  predict  a  model  whose  parameters  are  changing  rapidly  In  time  we 
may  require  a  more  responsive  system.  The  cost  of  this  Is  a  correspondingly 
higher  error  variance.  On  the  other  hand,  we  can  not  tolerate  an  infinite 
error  variance  because  this  Indicates  an  unstable  forecasting  system.  Such  a 
system  violates  exactly  those  assumptions  which  we  hold  most  Important  to  the 
generation  of  our  forecasts. 

3.2.  Practical  Consequences. 

It  Is  obvious  frdm  the  forms  of  V(T)  In  the  two  special  cases  (6)  and 
(7)  that  the  usual  range  of  values  of  smoothing  constants,  I.e.  (0,1),  lies 
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within  the  stability  regions  for  these  systems.  The  problem  arises  for  only 
the  general  seasonal  model.  The  actual  stability  region  1^  difficult  to  es¬ 
tablish  In  generel  here  for  It  depends  not  only  on  and 

on  n  .  the  length  of  the  season.  However.  Gardner  (1984)  reports  that  Sweet 
has  demonstrated  numerically  that  the  usual  range  for  the  three  smoothing 
constants  lies  within  the  stability  region  for  n  up  to  four.  For  seasons  ‘ 

longer  than  four  periods  It  Is  no  longer  true:  a  result  of  Importance  for 
weekly  and  monthly  data. 

The  values  of  the  smoothing  constants  are  established  In  one  of  two  ways 
In  general.  They  may  be  selected  Intuitively  by  appealing  to  the  Ideas  of 
required  speed  of  response  and  constancy  of  the  underlying  model.  Alterna¬ 
tively,  they  may  be  chosen  as  giving  the  dest  fit  with  respect  to  some  cri¬ 
terion  such  as  least  squares.  However  the  values  are  chosen,  we  shoqld  be 
very  concerned  If  they  do  not  lie  In  the  systems* s  stability  region. 

If  we  have  chosen  them  for  rapid  response  we  may  allow  their  use  for  a 
brief  period,  perhaps  when  we  Initialize  the  forecasting  system.  As  noted 
above,  however.  If  we  do  not  replace  them  with  values  In  the  stability  region 
then  the  longer  we  forecast  the  more  Influential  become  the  very  first 
observations.  If  some  best  fitting  criterion  selects  'unstable*  smoothing 
constants.  I.e.  ones  leading  to  en  Infinite  error  variance,  then  a  very  reel 
possibility  Is  that  the  model  Is  wrong,  at  least  for  part  of  the  data.  In 
particular.  It  may  well  be  that  a  multiplicative  seasonal  model  Is  called 
for. 

We  may  decide  to  retain  the  model  but  use  smoothing  constants  In  the 
stability  region.  These  will  be  obtained  by  decreasing  one  or  more  of  • 

Oq.  o^.  02  *  general.  If  o^,  o^,  02  do  npt  lie  In  the  stability  region 
It  will  be  because  one  or  more  of  them  is  too  large. 
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4.  SUMMARY 


The  purpose  of  this  paper  has  been  to  give  a  procedure  for  deriving  the 
variance  of  the  forecast  error  for  Winters'  Additive  Seasonal  System,  In 
summary  that  procedure  Is  as  follows: 

{1}  If  the  model  has  no  trend  or  Is  non>seasona1  the  results  are  given 
In  Section  2.1  and  2.2  respectively  In  dosed  form. 


(11)  From  the  chosen  smoothing  constants  Oq,  and  the 

coefficients  W^,  W^,...,  are  derived  as  In  Appendix  1. 

(Ill)  The  W^  values  are  used  to  start  the  algorithm  of  Appendix  2  and 
yield  dQ,d^,...,d^  .  The  algorithm  simultaneously  checks  the 
stability  of  the  forecasting  system  for  the  chosen  values 


oq.  02  •  If  n  _<  4  the  system  will  be  stable. 

(1v)  The  values  of  d^  are  substituted  Into  the  appropriate  expressions 
for  V(T)  and  Vg(T)  . 


Note  that  at  step  (1v)  we  are  able  to  evaluate  the  variance  of  the  errors  for 
T-step  ahead  forecasts  and  T-step  ahead  cumulative  forecasts.  Confidence 
Intervals  for  these  future  values  can  be  derived  from  these  variances. 

A  discussion  of  how  to  calculate  the  corresponding  results  when  the 
seasonal  factors  are  renoriMlIzed  after  each  observation  Is  given.  We  have 
also  discussed  the  (real)  possibility  that  the  error  variance  may  be  Infinite 
and  Indicated  Its  relationship  with  the  stability  of  the  forecasting  system. 
Some  recommendations. are  made  about  the  Interpretation  of  system  Instability. 
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APPEI^IX  1 


The  results  of  this  paper  are  obtained  using  the  fact  that  the  forecast¬ 
ing  system  given  by  (2)  has  an  equivalent  ARINA  process.  The  nature  of  this 
equivalence  Is  discussed  In  McKenzie  (1984)  and  the  equivalent  ARiriA  for  (2) 

Is  given  In  McKenzie  (1976),  It  has  the  form  (1-B)(1-b”)X^  »  W(B)ef.  ,  where 

k  k 

B  Is  the  backshift  operator  defined  by  B  X.  »  X^_|^  and  W(B)  ■  1  «  I  W.b'^ 

k*l 

where  *  1  -  Oq  ’  »o«i  *  “n  '  ^  “  ®0“l  ‘  *  “n+1  “  ‘  . 

and  M|j  ■  -  OyOi  ^  ®  2,3 . n-1  . 

Suppose  now  that  W(B)6(B)  *  1  where  6(B)  «  I  6.B*^  .  Using 

k-0  " 

equations  (2),  It  can  be  shown  that,  for  the  purposes  of  the  variance  calcula-r 
tion,  the  components  of  of  the  forecast  can  be  expressed  as  Infinite  moving 
averages  In  {a^}  .  Thus: 


"t  *  ^*0^^  ^  X  ^<^1  ■  » 

1“1 


m1n(k,n-l) 

where  ■  I  •  Further,  the  corresponding  representation  of 

1 ®U 

can  be  obtained  by  writing  It  as  B  S^  ,  (k  ■  1,2,,*, ^n)  . 

Since  {a^}  are  Independent  random  variables  the  variances  and  covariances 

can  be  obtained  directly  from  these  moving  average  representations.  Defining 

*  I  ”  U,l,...,n),  yields  the  expressions  given  In  the  paper. 

*  1»0 

The  algebra  Is  tedious  and  hardly  Illuminating  and  so  Is  omitted.  It  can  be 
obtained  upon  request  from  the  author. 
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As  rt^rds  tHt  rtfionis11<«t<oi;L  procedure »  recall  the  equivalent  ARIMA  for 
, ,  Xf  «!«,  Intro^fce  coatinuous  renormalization  as  discussed  we 

is,  fK|u1v«ltiit  APINA  and  it  has  the  same  form.  Now. 

1(0)  0^  replaced  by  aj,  oj,  .  This  ARIMA  Is  equlv- 

i1fflt>e  tiM.  syiteili  0)  with,  starred  smoothing  constants.  Hence,  the 
SMu)!. ,,  Is  surprising.  However,  note  that  we  are 

really  dealing  with  twd  different  decompositions  of  the  seasonal  factor.  One 
i,«fwil.ts  iS  It  ippfars  in  (2)  and  the  other  Is  where 

..  Is  normal  1  lid  t(f  sum  to  zero  and  m*  is  the  correspondingly  adjusted 

......  ,  ,  ^ 


’vr.  t  "i 


-  y.r^.r  -r’v'  *  ■'  ■  - 
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APPENDIX  2 

To  obtain  the  values  {d|^}  note  that,  by  definition,  they  are  the  variance 
(d^)  and  the  first  n  autocovarl ance$  aotoreyrefslve 

process  of  order  (n+1)  given  by  W(B)Z^  •  ,  ivhfere  fe^}  are  Independent 

random  variables  of  zero  mean  and  unit  variance*  Thus,  the  sequence  {d|^}  may 
be  obtained  by  solving  a  suitable  set  of  the  Yu1e>Wa1ker  equations  for  this 
process.  This  procedure  Is  discussed  In  McLeod  (1975,  1977). 

From  our  point  of  view,  however,  a  much  superior  approach  Is  presented  by 
Wilson  (1979).  The  stability  of  the  forecasting  system  corresponds  to  the 
statlonarlty  of  the  autoregressive  process  (Z^)  and  can  be  tested  routinely 
within  the  procedure.  The  algorithm  Is  as  follows: 

(1)  define  k  *  *  1,2,..., n+1;  and  t^^^j  »  1  . 

(11)  apply  the  following  equations  In  the  given  order  for 
k  =  n+l,n,...,2 

"k  •  *  -  "^.k 

If  D|^  <.  0  ,  the  system  Is  unstable;  stop  now. 

“k-1.1  *  <“k.)  *  “k.k  “k.k-i>^“k* «  ■  . '■-») 

‘k-i  •  ‘k/"k 

(til) 

If  0  ,  the  system  Is  unstable;  stop  now. 

dg  ■  t^/Dj,  and  the  system  Is  stable,  and  the  error  variance  finite, 

(t*)  'll  ■  >*1,1  <*0 


‘k  ■  “k-l.t  “k-l  *  “k.k  ‘k  •  ^  ■  2-’' 


Note  that  If  n  £4  the  system  Is  stable  for  the  usual  choice  of  parameters 
I.e.  c  (0,1),  1  *  0,1,2,  and  so  we  need  not  test;  D|^  . 


.  REFERENCES 


Bowerman.,  B,  L,  and  O'Connell,  R.  T,  Time  Series  and  Forecasting;  an  applied 
approach.  Massachusetts:  Duxbury  Press,  1^79, 

Box,  6.  E,  P,  and  Jenkins,  G,  M.  Time  Series  Analysis:  Forecasting  and 
Control  (2nd  ed.),  San  Francisco:  Holden-Day,  1976. 

Brenner,  J.  L.,  D'Esopo,  D.  A.  and  Fowler,  A.  G.  "Difference  Equations  in 
Forecasting  Formulas",  Management  Science.  15  (1968)  141-159, 

Brown,  R.  G.  Smoothing.  Forecasting  and  Prediction  of  Discrete  Time  Series. 
New  Jersey^  Prentice-Hall  1^62. 

Chatfield,  C.  "The  Holt-Winters  Forecasting  Procedure",  Appl,  Statistics. 

27  (1978)  264-279. 

Gardner,  E,  S.  "Exponential  Smoothing:  the  state  of  the  art",  J,  Forecasting, 
to  appear  (1984). 

McClain,  J.  0.  "Dynamics  of  Exponential  Smoothing  with  Trends  and  Seasonal 
Terms",  Management  Science.  20  (1974)  1300-1304. 

McKenzie,  Ed.  "A  Comparison  of  some  Standard  Seasonal  Forecasting  Systems", 

The  Statistician.  25  (1976)  3-14. 

McKenzie,  Ed.  "General  Exponential  Smoothing  and  the  Equivalent  ARMA  process", 
J.  Forecasting,  to  appear  (1984). 

McLeod,  A,  I,  "Derivation  of  the  Theoretical  Autocovariance  Function  of 
Autoregressive-Moving  Average  Time  Series",  Appl.  Statist*.  24  (1975) 
255-256,  (Correction,  Appl.  Statist.  26  (1977)  194.) 

Montgomery,  D.  C.  and  Johnson,  L.  A,  Forecasting  and  Time  Series  Analysis. 

New  York:  McGraw-Hill,  1976, 

Thomopoulos,  N,  T.  Applied  Forecasting  Methods,  New  Jersey:  Prentice-Hall, 

1980. 

Wilson,  6.  T.  "Some  Efficient  Computational  Procedures  for  High  Order  ARMA 
Models",  J.  Statist.  Comput.  Simul,  8  (1979)  301-309. 


17 


DISTRIBUTION  LIST 


Defense  Technical  Infornatlon  Center 
Caneron  Station 
Alexandria,  VA  22314 

Library 
Code  0142 

Naval  Postgraduate  School 
ifontOr^,  CA:  93943 

Research  Administration 
code  012A 

Naval  Postgraduate  School 
Nonterey,  CA  93943 

Library 
Code  55 

NafTal  Postgraduate  School 
Nonterey,  CA  93943 

Meteorology  Department 
Code  63 

Naval  Postgraduate  School 
TIooterey ,  CA  ^  93943 

Professor  Ed  McKenzie 
Code  55  • 

Navel  Postgr^ate  School. 

Monterey,  CA  93943 

Professor  P.'  1,  Leirfs 
Code  55Lw 

Naval  Postgraduate  School 
Morttereyj  CA  93943 


